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PLASMA EXTRACTION DEVICE 

This invention relates to the extraction of an ion beam 
from a plasma. The invention is particularly directed to 
5 the extraction of high brightness beams. Although there are 
several definitions of brightness, it is generally accepted 
that minimising emittance ensures the highest possible 
brightness for a given beam current. Emittance is a measure 
of the parallelism of the individual particle trajectories 

10 with- in a beam. In planar symmetry this means that 
particles follow rectilinear and parallel trajectories 
whilst in cylindrical and spherical symmetry they move along 
rectilinear paths that follow radial lines as though 
emmanating from or converging towards a single line or 

15 point. The present invention has application to diverging, 
converging and parallel ion beams extracted from rectangular 
slits and circular apertures. This corresponds to four 
basic beam types: 

Strip Beams: parallel beams extracted through a slit 
20 Wedge Beams: diverging/converging beams extracted 

through a slit 

Cylindrical Beams: parallel beams extracted through a 
circular aperture 

Conical Beams: diverging/converging beams extracted 
25 through a circular aperture. 

There are a range of applications for ion beams 
particularly in the semiconductor industry. For example, 
the fabrication and correction of lithography masks involves 
a need for sub-micron etching capability- This processing 
30 is currently achieved using medium-energy particle beams 
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(10-50keV) , most often ions, referred to as focused ion 
beams or FIBs. To ensure sub-micron feature creation, the 
FIBs must be f ocusable down to a nano-meter scale spot. 
This requires the extraction of very high brightness beams 
5 in excess of 10 5 A. sr^.m" 2 . Liquid metal ion source (LMIS) 
technology has been capable of this level of brightness for 
several decades. It exploits the capillary effect of liquid 
Gallium to cover a sharp Tungsten needle onto which a strong 
electric field is applied and consequently ions are removed 
10 by field effect. The field effect is strongest at the 
needle point and so a beam of ions is created that appears 
to emanate from a nanometer spot. This beam is then 
accelerated and focused onto the target where it etches the 
surface by surface collision processes. 

Though this technology boasts a nanometer scale milling 
capability, it produces nefaterious doping effects by 
introducing Gallium into Silicon. To avoid this a high 
brightness beam of inert gas ions would be preferable. This 
could be achieved by extracting ions from an inert gas 
plasma. However^ this has proven notoriously difficult and 
has been the subject of intense- study for some decades. 
Much work has been done to optimise extraction electrodes in 
a number of ways . For example adjustment of aperture ratios 
and electrode spacings is described in J.R. Coupland, T.S. 
Green, D.P. Hammond, and A.G. Riviere, A study of the ion 
beam intensity and divergence obtained from a single 
aperture three electrode extraction system, l . Rev. Sci . 
Instrum. , 44(9):1258, 1973. Shaping of electrodes has also 
been investigated as described in D.E. Radley, The theory of 
the pierce type electron gun, J. Electron . Control , 4:125, 
1957, E.R. Harrison, Approximate electrode shapes for a 
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cylindrical electron beam, Brit. J. Appl . Phys . , 5:40, 1953 
and P.N. Daykin, Electrode shapes for a cylindrical electron 
beam, Brit. J. Appl. Phys., 6:248/ 1955. 

Despite these efforts brightness in excess of 10 5 A. sr* 1 . m" 2 
5 has not been achieved. 

Accordingly this invention provides for a plasma extraction 
device for extraction of a high brightness ion beam of 
selected profile said device including at least three 
sequentially arranged electrodes each respectively having an 

10 aperture corresponding to the selected beam profile, a first 
of said electrodes adjacent said source being maintained at 
an extraction potential, a second of said electrodes being 
spaced from the first electrode and maintained at a first 
blocking potential opposite to said extraction potential; 

15 and a third of said electrodes disposed between said first 
electrode and said second electrode and maintained at a 
selected potential between said extraction potential and 
said first blocking potential; said first, second and third 
electrodes being shaped to produce substantially zero 

20 electric field at said third electrode. 

Preferably, the extraction device includes a further two 
sequentially arranged fourth and fifth electrodes each 
respectively having an aperture corresponding to the 
selected beam profile and disposed on the side of said 

25 second electrode remote from the plasma, said fourth 
electrode being between said second and fifth electrodes and 
maintained at a second blocking potential less than said 
first blocking potential and said fifth electrode being 
maintained substantially at zero potential, the fourth and 

30 fifth electrodes being shaped to maintain substantially zero 
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electric field at the second and fifth electrodes. 

In preferred forms of the invention the beam profile is one 
of a strip beam, a wedge beam, a cylindrical beam, and a 
conical beam, 

5 These beams can be thought of as sections of current flowing 
in complete diodes. Consequently, parallel beams can be 
treated in the same way as a planar diode, 
divergent /convergent beams extracted through a rectangular 
slit can be treated as a cylindrical diode and 

10 divergent /convergent beams extracted through a circular 
aperture can be treated as a spherical diode. In each case, 
the plasma/beam interface or plasma meniscus is the anode 
(emitter) and the 0V equipotential surface is the cathode 
(collector) . The purpose of the extraction device is to 

15 ensure that the meniscus and 0V equipotential surfaces are 
parallel in the case of parallel beams and concentric 
cylinders or spheres in the case of diverging or converging 
beams. If this situation is maintained then, neglecting the 
inherent ion temperature in the plasma, the ions travel in 

20 perfect parallelism along radial lines from the meniscus to 
the 0V equipotential and suffer no deflection. In other 
words the beams have zero emittance growth. 

In accelerating the beam, however, a significant electric 
field can be produced in the direction of flow. In 
25 accordance with the present invention this is compensated to 
avoid the electric field in the beam from ballooning 
outwards at the exit of the extractor and deflecting the ion 
trajectories. 

In the transport region aft of the extractor, considerable 
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advantage can be gained by the presence of electrons which 
neutralise the beam space charge* However, their presence 
in the extractor, especially in th^ acceleration gap, is 
highly detrimental to beam quality since they alter the 
5 charge density distribution and hence the potential 
structure through-out the extractor. The existence of an 
electric field at the exit to the extractor would serve to 
accelerate electrons into the extraction region in such 
quantity as to neutralise the electric field at the 
10 extractor exit. Accordingly the present invention ensures 
that the beam potential has zero electric field upon. exit of 
the extractor. 

The electrons in the beam plasma aft of the extractor have 
non-zero temperature and are generally distributed according 

15 to Maxwell's law* This means that higher energy electrons 
from the tail of the distribution are able to enter the 
extractor if a blocking field is not present. To this end a 
small blocking potential preferably of some hundred Volts is 
produced at the end of the extraction region to inhibit the 

20- passage of. electrons. Again, to avoid ballooning of the 
electric field this potential must be produced so that the 
electric field upon exit of the blocking region is zero. 

In a preferred form of the invention the potential 
25 distribution in the beam will then have the form shown in 
Figure 1. There are two regions called the Extraction and 
Blocking regions both of which are divided into two stages. 
Stage 1 serves to take the beam from a low gradient to a 
high gradient and stage 2 vice versa. In the Extraction 
30 gap, stage 1 is necessary to match the beam potential to the 
plasma sheath and stage 2 to bring the electric field to 
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zero. A similar rationale is behind the two stages of the 
blocking region. 

The extractor region is formed by the first three sequential 
electrodes. Stage 1 and stage 2 are separated by the third 
5 electrode. The first electrode is preferably the plasma 
electrode. The second electrode and third electrode are 
usually referred to as the blocking and accelerating (Accel) 
electrode respectively. 

The blocking region is formed between the second electrode 
10 and the fifth electrode. Stage 1 and stage 2 are separated 
by the fourth electrode. 

In order to describe the plasma extractor device of the 
present invention it is firstly necessary to provide some 
analysis of beam distributions and an overview of solutions 
15 to Laplaces equation for different beam profiles. This will 
be done in separately headed sections which proceed an 
explanation of electrode design and specific examples of the 
invention. Reference will be made to the accompanying 
drawings in which: 

20 Figure 1 is a schematic graph of potential distribution in a 
beam extracted according to the present invention; 

Figure 2 is a schematic representation of a diverging beam 
showing the extractor and blocking stages according to the 
present invent ion ; 

25 Figure 3 is a graph of generalised distribution versus 
Langmuir-Blodgett ; 

Figure 4 illustrates a contour C used in equation (46) ; 

Figure 5 is a plot of electrode geometries forming a plasma 
extraction device according to the invention for a 5° 
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diverging wedge beam; and 

Figure 6 is an enlargement of part of the plot of Figure 5. 

Figure 1 shows the general shape of the potential along the 
beam inside the extractor for low emittance growth 
5 extraction (zero emittance growth in the ideal case) . The 
details of these distributions in each stage of the 
extractor are first explained. 

The overall potential profile shown in Figure 1 is in fact 
four distributions concatenated. Each stage is bounded by 

10 two plane or two concentric surfaces depending on the beam 
type, so that each stage can be treated as separate, 
•complete' diodes, each with its own set of boundary 
conditions. To get the complete distribution, the solutions 
are stitched together by matching the boundary' conditions at 

15 each surface. 

The potential and charge distributions in the beam at each 
stage of extraction are governed by Poisson's law. The 
basic solution to this problem, using simple boundary 
conditions, in plane symmetry was solved by Child and 
20 Langmuir in 1911 and 1914 respectively and in cylindrical 
and spherical symmetry by Langmuir and Blodgett in the 
1920s. These initial solutions assumed only one charged 
particle species (notably electrons), and ignored initial 
velocity. 

25 In accordance with the present invention several 
generalisations are introduced. In particular, 

distributions describing non-zero initial gradient, 
distributions describing non-negligible initial velocity, 
distributions tapering from a strong gradient to a zero 

30 gradient (the reverse of the standard Langmuir-Blodgett 
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solution) and distributions describing the presence of 
Maxwellian electrons. A further requirement on these 
distributions is that they have the same form as the 
standard Langmuir-Blodgett series solutions since the Radley 
5 solutions to Laplace's equation are dependent on this form. 

The work published by Langmuir-Blodgett in the 1920s forms 
the basis for the beam distribution analysis that follows 
and is given in terms of simple and compact series 
solutions. A cursory mathematical overview of Langmuir- 
10 Blodgett 1 s contribution starting with spherical symmetry and 
working through cylindrical to planar symmetry is provided. 

In the case of ion beam extraction from a circular aperture, 
the charge and potential distributions in the beam are 
assumed to be analogous to those in a complete spherical 
15 diode. Following Langmuir and Blodgett, Poissons 1 equation 
between two concentric spheres can be stated as 




where V is the potential at a point distant r from the 

common center and p is the ion charge density. The current 

20 flowing in the diode can be written in terms of the particle 
velocity v m > 

/=Wpv (2) 

where the velocity can be written in terms of the voltage V 
by using the kinetic energy relation: 

25 -Mv 2 =-eV (3) 

2 

Combining equations (!) , (2) and (3) yields: 



m 
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% d 1 V , rfK , f7 -i 

r — -r + 2r ~A V 2 (4) 

dr dr 

where V =-V , and: 

A.J-M (5, 

Equation (4) can probably not be integrated directly but a 
5 solution can be found in terms of a series. The form of the 
solution as a function of the ratio K=r/r 8 



V(R)=\^A\ /»■(*) (6) 

(9 V 

where / is the analytic function to be found. The -Ir^J 

term serves to normalise for the constant term A related to 
10 the current and hence the plasma density and meniscus 

4 

curvature, and the f*(R) term to remove the square root from 
($) and hence to simplify subsequent derivations. One 
further transformation is performed by setting: 



r = ln(i?) (7) 



15 so that a solution to (6) can be expressed in terms of a 
MacLauren series : 



f = £a„y n (8) 



substituting (6) into (4) and imposing equation (7) results 
in: 

20 + +3^-1 = 0 (9) 

where f'-dfldy and f" = d 2 f/dy 2 . From this form it can be 
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seen that for / = 0, /' = !- The first six terms of the series 
solution are then: 

/«^-0.3^ 2 +0.075^ 3 -0.0143182r 4 +0.0021 609r 5 -0.0002679 If 6 (10) 
where this corresponds to the case K~0 and K' = 0. 

5 A similar derivation can be made for the case of cylindrical 
symmetry. Poisson's equation becomes: 

dr dr 

where : 

g= 1 r— (12) 



2* 

10 and 1 is the length of the extraction slit- The solution 
takes the form: 



■(13) 



where g is the analytic function to be found oyer the 

desired range of r . Again substituting (13) into (11) and 
15 setting (7), results in: 

3£g" + g' 2 +4gg' + g 2 -l = 0 (14) 

Again, for g=0 f = 1 and the series becomes: 

g = r~0Ay 2 +0.0916667y 3 -0.01424242r 4 +0.00 1679275^ 5 -0.000 16 122 19^ 6 (15) 
That this corresponds to the case V = 0 and K' = 0. 
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2.2 Non-Zero Initial Gradient 

The present invention is direct to the extraction of ions 
from plasma, the assumption that the electric field at the 
meniscus surface is zero, central to the Child-Langmuir and 
5 Langmuir-Blodgett derivations, is incorrect and so will be 
generalised. Furthermore, the final result will be 

presented in a similar, easy to use series formulation as 
the original Langmuir-Blodgett relation. 

Figure 3 demonstrates the effect of non-zero initial 
10 gradient. The dashed line corresponds to the case of a 
divergent conical beam extracted from a plasma of density 
K = 10 M cm" y and is compared to the standard Langmuir-Blodgett 
distribution which assumes zero initial gradient. 

To understand why the gradient at the plasma/sheath 

15 interface is meniscus is non-zero it is necessary to 
consider the Bohm sheath criterion which stipulates the 
minimum ion velocity for entry into the sheath to maintain a 
stable sheath at a plasma boundary. In conjunction with 
. some distribution relation for electrons, this defines a 

20 potential structure with-in the sheath. In particular, the 
electric field at the plasma boundary is non-zero and is 
typically several hundred kilo- Volts per meter. For 
continuity of the electric field across the plasma/beam 
interface the potential gradient must be equal on both sides 

25 of this surface. This is not assured by the original 
assumptions of Langmuir and Blodgett who were modelling 
particle flow from thermionic cathodes. In that case the 
source of particles was assumed to be undeletable and to 
have no intrinsic electric field so that for equilibrium the 

30 boundary condition in the extractor was for the electric 
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field to be zero at the entry to the acceleration gap. In 
the case of plasmas , the flux of ions is fixed and an 
intrinsic" electric field does exist in the sheath which 
separates the beam from the bulk plasma. This strongly 
5 implies that particle beam extraction from plasmas is not 
space charge limited but rather is source limited and 
further that the voltage distribution in the beam is not 
given by the Langmuir-Blodgett relation as it is usually 
stated. 

10 General solutions to Poisson's equation for the case of non- 
zero initial gradient are first present followed by a 
discussion of the plasma sheath and how the gradient at the 
meniscus is obtained. 

Spherical Symmetry 

15 Taking the first derivative of (6) , yields: 

r{y)^x4p (16) 

1 

where A=^^/lj 3 and p = f'{y) 2 fiy)* * for y = Q. In the classic 

Langmuir-Blodgett derivation /(^) = 0 in (9) leads to /'(^) = 0, 
for y = \ f (assuming potential increases as a function of 

20 position in the extractor), so that p = 0 and hence 
V'{y = Q) = 0. In other words, the Langmuir-Blodgett derivation 
requires a zero initial gradient in potential. Numerically, 
however, it is possible to have f approach zero without 
requiring K'(o) to annul by setting /'(o) such that (16) holds 

25 for the desired value of V'(®) • For a given value of K'(o) and 
A , there is a limit to how small /(o) can be set, but in 
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most practical cases it is several orders of magnitude less 
than unity. 

The solution depends only on p, rather than the individual 
values of f(0) and f'(0) , for the range of p which is of 
5 interest. The series coefficients, a n ,are expressed as a 
quadratic: 

a H =a n +P n p + YnP 2 (17) 

where a n , fi„ and y n are the expansion coefficients found by 
a least squares method. These terms are presented in the 
10 following table 1. 



n 




A 


r„ 


1 


1.0035 


4'. 0.49 


-10.92 


2 


-0.3084 


-8.008 


25.11 


3 


0.08338 


7.791 


-25.85 


4 


-0 . 01825 


- 3 . 96 


13.47 


5 


0. 002870 


1 . 0004 


-3 . 448 


6 


-0. 0002227 


-0.09904 


0.3441 



Table 1: Expansion terms for the coefficients of the 
MacLauren series (see equation 17) . Spherical case. Note: 
15 the a n are very close to the original Langmuir-Blodgett 
series coefficients. 

Cylindrical Symmetry 

Taking the first derivative of the potential and writing it 
in terms of the parameter p : 
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(18) 



where ft 



and /J=l, and note that in the limit the 



-KM 

term in g disappears, so that V'(R)*>2jjJp This has the 
same form as the spherical case, so that by writing: 



(19) 



and plotting the series coefficients in terms of the 
parameter p yields: 



(20) 



where 5 n , £*„ and are coefficients found in the same 

10 fashion as for the spherical case. Results are presented in 
table 2 . 



n 


Yn 


8. 


e » 


1 


1 . 0034 


3 . 989 


-10,69 


2 


-0 .4086 


-8.223 


25 .43 


3 


0 . 1005 


7.974 


-26.22 


4 


-0,01866 


-4.038 


13 .65 


5 


0.002658 


1.0184 


-3.491 


6 


-0 . 0002011 


-0. 1007 


0.3482 



Table 2 : Expansion terms for the coefficients of the 
MacLauren series (see equation 20) ♦ Cylindrical case. 
15 Note: the y n are very close to the original Langmuir-Blodgett 
series coefficients . 
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The Plasma Sheath 

There are several ways to model the sheath. . A method in 
terms of the Bohm sheath criterion and the Boltzman equation 
in one dimension, which is easily extended to spherical and 
5 cylindrical symmetry is described in full in M. A. Liberman 
and A. J. Lichtenberg "Principles of Plasma Discharges and 
Materials Processing" John Wiley & Sons, New York, 1 st 
edition, 1994. This model incorporates both the non-zero 
ion velocity v B at the entry to the sheath, required by the 

10 Bohm criterion, and also the presence of electrons. Another 
popular method is the Child sheath, which is extended to 
spherical and cylindrical symmetry by the use of the 
standard Langmuir-Blodgett corrections. In this case, the 
pre-sheath/sheath boundary and the meniscus are considered 

15 to be concentric spheres for extraction from a circular 
aperture and concentric cylinders for extraction from a 
slit. Though solving the Boltzman sheath is possible in 
terms of a series, it requires a somewhat more drawn out 
analysis and so for simplicity, the Child sheath method will 

20 be employed with Langmuir-Blodgett corrections. As such, 
the pre -sheath can be ignored and we assume that the 
velocity of ions and the potential at the bulk plasma/sheath 
edge are zero. However, it should be noted that the Child 
sheath yields smaller gradients than the Boltzman sheath. 

25 Sheath Potential at the Meniscus 

To solve (4), the case of extraction fro a plasma requires / 
must equal the ambipolar flux for ions: 



/* =0.6en i - v a - A 



(21) 
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where n s is the plasma density at the sheath edge and A = 4^r, 2 
is the area over which current is extracted, r, being the 
radius of curvature of the sheath. It follows from (5) and 
(21) , that the solution to (4) is strongly relates to both n t 

5 and r s 

In addition, three boundary conditions are required.* At the 

entry to the sheath set P(l)=0 and = 0 ■ To determine the 

sheath potential at the meniscus equate ion flux, assumed 
constant through-out the sheath 

io r ' = ^ (22> 

to the electron flux at the meniscus, 

r « 4/? * (23) 

where = {8eT e 1 7un)1 is the mean electron velocity and V m is 
the potential of the meniscus with respect to the plasma- 
15 sheath edge. Thus upon substitution of the Bohm velocity: 



20 



which becomes: 

i 



This can be expressed in a more convenient form by 
substituting the mass of the extracted ion species. Krypton 
is a typical gas used in ion beam extraction from plasmas 
and has a mass of M - 84<w . Therefore in this case, (25) can 
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be rewritten as V m = -5,052T e , which, assuming T e = 3eK , is 
approximately -15V. This now leads to suitable boundary 
conditions for (4) : 



V = -\W ,r = r m (26) 
dVidR = 0 >r = r, 



where r„ is the radius of curvature of the meniscus. 
Potential Gradient at the Meniscus 

Since both the voltage and its first derivative are zero at 
10 r gt f is independent of the parameter p, from the familiar 
Langmuir-Blodgett relation: 

a( r )=^-0.3^ 2 +0.075r 3 -0.0i43l82/ 4 +0.0021609r 5 -0.00026791/ 6 (27) 

15 This series expansion in conjunction with (6) and the 
boundary conditions now determine both the sheath width and 
the potential gradient at the meniscus edge. The sheath 
width is taken as the value of r^-r^ for which (6) is equal 
to (25). The potential gradient at the meniscus edge is 

20 then equal to the first derivative of (6) taken at this 
value of r m . 

Assuming a constant electron temperature and gas type, 
equations (6), (7) and (27) show that the potential gradient 
at the meniscus edge is dependent on the bulk plasma density 
25 and the radius of curvature of the meniscus . 



Solving Poisson's Law Backwards 
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To solve this consider the standard Langmuir-Blodgett 
problem in reverse, solving the differential equations from 
the exit of the extractor where the gradient is zero to trhe 
entry where it is large. This is tantamount to reversing 
5 the distribution in the case of a parallel beam or solving 
for the opposite convergence in the case of divergent or 
convergent be ants . 

The f and. gr series as defined by 10 and 19 remain unchanged 
however the definition of J* and hence y is altered. In the 

10 case of a diverging beam, R is taken to be the ratio of the 
current position to the first concentric surface and thus 
greater than unity. But* in the case where Poisson's 
equation is solved backwards it is redefined to be the ratio 
of the current position to the second concentric surface and 

15 hence less than one, The inverse is true for a convergent 
beam. 

Presence of Electrons 

To account for the presence of a population of electrons 
arriving from the tail of a Maxwellian distribution an 
20 exponential term is added to the original differential 
equations to account for the Boltzman relation ♦ . 

Generalising for the presence of electrons involves only the 
RHS of Poisson's law: 

p e(n, -n e ) 

— = — — (28) 

25 Since the electrons belong to a Maxwellian distribution 
their density as a function of potential is given by 
Boltzman' s law: 
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(29) 



where n 0 is some percentage ^ of the ion density in the beam 
aft of the extractor to take into account the fact that 
neutralisation is not always 100%. The differential 

5 equations for spherical, cylindrical and planar symmetry 
then become : 



^^ + Cexp(C)j f 



r 2 K* + 2rr = 0.6wd^ + ^exp(C) , where A=v B A, 



2e 



(30) 



rV* + V = 0.6ne£jL + ^ exp(C) 



where B = v» J— 
* V2e 



10 when n is the ion density in the beam, A is the area of the 

eV 

anode and C= . Here the solutions take a very different 

kT e 

form to those presented previously on account of the 
exponential term. However, it is possible to give a 
solution of the potential in terms of a MacLauren series, 
15 which when suitably normalised yields the correct form for 
implementation in the solutions to Laplace's equation. 

Non-Zero Initial Velocity 

The problem is set out in the same fashion as in section 2.1 
for negligible initial velocity, except that now the kinetic 
20 energy relation is written: 

Lmv 2 ~~Mvl=eV (31) 
2 2 
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This when rearranged yields: 



V + v 2 0 (32) 



But since: 



v2 °=$F"' (33) 

where V al is the energy of the particles at the exit of the 

previous stage, the differential equations stipulated by 
Poisson's law become for spherical, cylindrical and planar 
symmetries : 

r 2 V" + 2rV'=- A 



10 f + K « (34) 

And the solutions to these equations become: 



(9 Y> - 

w 3 - 



(35) 

~ 4 



-(Hv-"- 



Importantly, upon substituting (35) into (34), (9) and (14) 
remain unchanged. This means that the initial velocity 
15 serves only to translate the solution vertically. The f and 
g series can still be obtained by the various means set out 
previously. 
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Laplace v s Equation 

In a given region of the extractor according to this 
invention/ determining the electrode geometry amounts to 
solving Laplace's equation subject to the potential along 
5 the beam edge. Since in the three cases of interest, plane, 
cylindrical and spherical geometry there is strong symmetry, 
this means that compact solutions can be obtained. These 
were presented by Radley [3] in 1957 along with a complete 
and rigorous mathematical derivation. A cursory overview is 
10 provided in the following. 

However, before doing so it is noted that in treating the 
instability issues of the solution, Radley remarks that in 
as much as small variations of the initial surface can 
produce large differences in the solution so, conversely, do 
15 relatively large variations in electrode shapes away from 
the beam surface produce only small variations in the form 
of the beam surface. This will have important implications 
in the section Electrode Design because electrodes will have 
to be curtailed to avoid overlap or break-down proximity, 

20 Strip Beam 

A strip beam can be thought of as an infinite plane diode in 
which the cathode is the plane x=0 and all charge in the 
region y>0 has been suppressed. To determine electrodes 
that would extract a beam of this sort a family of 
25 equipotentials in y>0 must be found such that the conditions 
in y<0 are unchanged. The basic case of the potential 
distribution in space-charge limited flow when suitably 
normalised yields: 
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4 

K = ** , in y^O (36) 

Since, assuming zero ion temperature, the trajectories are 
rectilinear and perpendicular to the emission surface 
defined by the extraction slit, the transverse forces on the 
beam at the beam edge are zero, so that 5Vldy = Q there. 
Restating the conditions in polar coordinates, gives: 



4 

^ = 0 
89 



on 0 = 0 (37) 



The solution in Cartesian coordinates is: 

K(jc,y)«Rcj|z. 7 | (38) 



10 where 

z = x + iy-r>e 10 (3 9) 

so that the solution is given by the conformal mapping 



4 

W = V + iU = (40) 



Wedge Beam 



15 For a wedge beam it is assumed that the meniscus and 0V 
equipotential surfaces are concentric cylinders. The 
coordinate system is therefore chosen so that the origin is 
at the vertex of the wedge, and the beam surface lies on 
0-0. Since the trajectories are rectilinear, lying along 

20 the lines 0=const., SVISB^O on 0-6. The meniscus is taken 
to be r=r m and the cathode to lie somewhere in the region 
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r>r M in the case of divergent beams and r<r m in the case of 
convergent beams* In the following it is assumed that the 

ratio /?=— is greater than unity or that the wedge is 
r m 

divergent. The potential distribution along the beam edge, 
5 when suitably normalised, is given by: 

2 4 

V = R*fP (41) 
where p is an infinite power series in the variable y = ln(/?) . 
(41) can be rewritten in series form: 



(42) 



10 Since a t = 1 , this last expression may be expanded by the 
multinomial theorem, to give: 



2 A 



v = R y y*Y, b «r n ~ ] (43) 



The coefficients Jb rt will be presented and discussed later. 
. The potential outside the beam is obtained by writing Re w 
15 for R in (41) , Thus y is replaced by co = y + i0, so that : v 



2 l i9 1 « 1 
R" <e> Q)^b n co n ' 1 (44) 

H«l J 



K = Re 



Cylindrical Beam 

For this case cylindrical polar coordinates (r,6,z), are 
enlarged with the axis of the beam along r=0. By suitably 
20 normalising the coordinates, the beam surface can be taken 
as r=l. Since the system has axial symmetry, none of the 
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variables involves 0, so this coordinate may be neglected. 
The ion flow considered is a cylindrical section of an 
infinite plane diode so that: 



on r=l (45) 



5 Applying the analysis in Radley, the solution is: 

7V 1 , 4 . s a ^^^w-^v.wvf <«> 

where p is a complex parameter, J and Y are Bessel functions 
of the first and second kind and C is the contour defined as 
shown in Figure 4. 

10 In reality, the contribution of both straight line segments 
in C cancel and the contour reduces to a circle of radius p . 
In theory any value of p will work but it has been found 
that values of between 1 and 5 produced the most rapid and 
accurate results. 

15 Conical Beam 

Spherical polar coordinates (r, 0,$)-, are employed with a 
cone semi -angle of 0=0 O , the anode on the sphere <R=1 and the 
cathode in the region J?>1 in the case of a diverging beam 
i and R<1 in the case of a converging beam. Again a diverging 
20 beam for this is assumed derivation. Writing y = ln(/?) and 
suitably normalising the potential gives: 



4 

v = p 



(47) 
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where 



/ = £=c.r" - (48) 



The coefficients c n depend on the boundary conditions for the 
beam and are discussed elsewhere. Again this series 
5 representation is expanded using the multinomial theorem and 
becomes : 



n-l 

The boundary conditions to be applied are: 



on 0=0 O (50) 



10 and following Radley gives: 

■•'K-^-'I-i) - 

(51) 

where u is a complex variable, P 0 and Q tf ^ re Legendre 
functions of the first and second kind and C is the contour 
.15 defined in Figure 10. Again the straight line segments 
cancel so that the contour reduces to a circle, however, now 
the radius p is dependent on R* 

The first step in the design of electrodes is the choice of 
beam type, beam current and final extraction energy. For 
20 each choice there is a different electrode design. It must 
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be determined from the first instance whether the beam is to 
be strip/wedge or cylindrical/conical and whether it is to 
be parallel, convergent or divergent. It must be kept in 
mind that the whole extractor is to act as if the beam was 
5 part of an entire diode from the plasma to the OV 
equipotential . The electrode design is such that the 
extractor apertures follow the beam shape so that they are 
just in contact with the beam at their respective positions. 

The Extraction Gaps 

10 Stage 1 

The extraction gap stage 1 comprises the aft face of the 
Plasma electrode and front face of the Accel electrode. 
Determining the beam distribution in this region, first 
takes into account that the electric field at the 

15 plasma/beam boundary is non-zero due to the plasma sheath, 
then given the beam form required the relevant solution can 
be arrived at by following the derivations above. This in 
turn is substituted into the relevant solution of Laplace's 
equation. Note that in this region the potential 

20 distribution is convex and that the electric field at the 
exit of this region is very large. 

Stage 2 

The extraction gap stage 2 comprises the aft face of the 
Accel electrode and the front face of the first Blocking 
25 electrode. The purpose of stage 2 is primarily to bring the 
electric field at the exit of stage 1 to zero. This is 
necessary as failure to do so will result in a strong 
ballooning of the electric field aft of the Accel aperture. 
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This will induce strong aberration and readjustment of the 
beam charge distribution leading to non- uniformity in the 
beam and potentially also reshaping of the meniscus away 
from the ideal plane, cylinder or sphere. The potential 
5 distribution in this region is obtained by solving Poisson's 
law backwards and assuming a non-negligible initial 
velocity. Again this expression is substituted into the 
relevant solution of Laplace's equation, 

Laplace's equation, it seems, can not be solved (in the Real 
10 domain) for a boundary condition that changes convexity 
because this would require that the equipotentials overlap, 
which in terms of electrodes means that they would need to 
occupy the same space. In fact, this is only a major 
problem for the Accel electrode and a solution is to shape 
15 the electrode so that it is the median between the two ideal 
cases. It should be noted that close to the beam, the two 
equipotentials are almost identical and that away from the 
beam they are not strongly disparate. 

Though the outer electrodes would eventually also overlap, 
20 and in practical terms would approach each other so that the 
inter-electrode gap would lead to break down, Radley has 
indicated that the effect of the electrodes away from the 
beam edge is increasingly negligible. Thus as a best 
approximation to an ideal extractor, the outer electrodes 
25 are made to extend to just outside break-down distance and 
the intermediate electrode is made to be the average of the 
two ideal equipotentials. 

Blocking Electrodes 

The Extraction electrodes have accelerated the beam just 
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beyond the desired extraction energy and have done so in 
such a way that the electric field at the exit of stage 2 is 
0 Volts per meter: In the case were positive ions are being 
extracted, this means that the voltage .at the exit of the 
5 Extraction region stage 2 is some negative value. The 
purpose of the Blocking electrodes is then to bring the beam 
potential back up to OV whilst ensuring that the electric 
field upon exit is 0 Volts per meter. Assuming that this 
Blocking potential is sufficiently high this impedes the 
10 flow of electrons from aft of the extractor to the 
Extraction region. 

Since the potential gradient is zero upon entry to the 
Blocking electrode stage 1, the standard Langmuir-Blodgett 
representation can be used for the beam potential. For 

15 stage 2, however, the presence of electrons can not be 
ignored since a population of higher energy electrons from 
the tail of the distribution will be able to penetrate some 
distance up the potential well. To solve for this 
distribution the relation given above is employed. This 

20 potential distribution is then substituted into the relevant 
solution of Laplace's equation. 

Beam Neutralisation 

Beam neutralisation is necessary aft of the extractor to 
compensate the considerable beam space charge. In many 

25 practical solutions, secondary electron emission from 
sputtering will be sufficient to provide a population of 
neutralising electrons. However, if this were not the case, 
some electron source would be required such as a hollow 
cathode. In this situation emphasis should be given to 

30 ensuring that the electrons are as low energy as possible to 
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avoid needing a large Blocking potential. 

. Plasma Density Range _ 

There is a limit to how high the extraction energy can be 
taken as a function of density and beam form. For example 
5 at a plasma density of 10 14 cnf 3 , no wedge beam form dan be 
extracted without aberration because the minimum voltage per 
meter required for extraction is in excess of the breakdown 
limit of 10 7 V/m. At 10 13 cnf 3 extraction of all beam forms is 
possible, but for convergent beams the maximum extraction 
10 energy is capped at 20kv in the stage 1 of the extraction 
region. Above - this energy, the electrodes need to be too 
close for breakdown. 

Examples 

Diverging Wedge Beam 

15 For this example a 5° diverging beam is extracted through a 
1mm wide, lm long rectangular slit. The plasma source is 
assumed to be a Krypton plasma of density 10 13 cm° at the 
sheath edge. The total extraction energy is 20kV and the 
total blocking potential 200V. 

20 Extraction Gap Stage 1 

The problem can be considered analogous to that of a 
complete cylindrical diode of curvature r*=5. 737mm. The 
current density flowing across the meniscus yields 
J3=7 . 6218- 10 8 according to (12). The ratio of the sheath to 
25 the pre- sheath is given by first solving (6) for the 
boundary conditions set in (25) : 
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P{y„)= , l5A ■ = 1.71827 -10- 3 (52) 



and then solving (19) to find y m =l. 71797 -10* 3 . From (7) 
K=l. 00172 so that the sheath width is given by: 

r, 1.001 72 -r, =9.9pm (53) 
5 The gradient at the sheath edge in terms of y m is: 

~{y m = 1.71827 10- 3 )=11644^/«m7 (54) 
ay 



Since y = lnf — I, dy=— Which 



means : 



^1 ( r > 1 " « 2.03 ■ 1 0 V.*.- (55) 

According to (18), p=l . 43557 • 10" 2 which is well below the 
10 established limit of Piimit=0.15, In combining table 2 and 
equation (20) the series expansion of gr is found to be: 

1.05846r-0.521406x 2 + 0.209569r 3 -0.0738154/ + 0.0 1 65584^ 5 -0.001 57496^ 6 (56) 

and hence an expression for the voltage in terms of y 
through relation (13) . Using the multinomial theorem (56) 
15 can be rewritten to give the form: 

1.0787 -0.70850 1^ + 0.342937;' 2 -0.140694^ 3 -kO.041 9959/ 4 - 0.00842 186? 5 (57) 



Since the wedge is divergent y = ln(/?) , and we solve (44) with 
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*=.»-■ ' ■ (58) 



where the non-normalised' coordinates are ; 



= r m 3c 



(59) 



For the aft face of the plasma electrode we solve (44) for 
V=0. To solve for the front face of the Accel electrode we 
solve for: 

V= l0AkV =0,138575 (60) 



(W 



which is a direct consequence of (41) . 
Extraction Gap Stage 2 

At the entry to stage 2, V ext =10 • lkV . The solution to the g 
series is given by the standard Langmuir-Blodgett relation 
but derived backwards ; Hence substitute: 

g « Y + 0.4/ 2 +0.0916667x 3 + 0.0 1424242/ 4 +0.00 1679275/ 5 +0.00016I2219/ 6 (61) 
into equation (44), but with ^ = -ln^j where r d is the radius 



of the concentric surface at the exit of the stage. At this 
point in the calculation is unknown, but by ( an iterative 
process a value can be found such that the potential 
gradient at r a at the entry to the stage is equ^l to that of 
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the distribution at the exit of the previous stage. We 
define: 



Ya =" ,n 



(62) 



so that 



(63) 



and: 



(64) 



Solve for y a so that 



= 63586 



(65) 



10 noting that the gr series in stage 1 is different to that in 
stage 2. Thus y 0 =0.35464. The voltage at y G is V=. 239885 
according to (41) . 

The Decel 1 electrode front face is calculated by solving 
(44) for V=0 and the Accel aft face is calculated by solving 

15 (44) for V=. 239885. It should be noted that the 

acceleration voltage in this gap is 174 84V which is 
significantly higher than lOlkV. To achieve a total 
acceleration energy of exactly 20kv, an iterative approach 
will be required. But since this does not benefit the 

20 illustration of the method, this will not be done here. 



Because the particles have an initial energy of lO.lkV, the 
basic solution is shifted up by this amount. 
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Blocking Gap Stage 1 

A blocking voltage of 2 00V is required so it is .assumed that 

the voltage in this stage is 10 0V. The standard Langmuir- 

Blodgett relation is solved to determine the gap size. This 

5 yield y =.00712381. Then to solve for this stage replace the 

standard Langmuir-Blodgett relation into (44) , For the 

Decel 1 aft face V=0 is solved and for the decel 2 front 

i 

face V=. 00137203 is solved. The final solution is shifted 
up by 27584V. 

10 Blocking Gap stage 2 

The solutions for the Poisson equation with the presence of 
Maxwellian electrons has not been resolved. For this reason 
this analysis does not assume the presence of electrons* 

The gradient at the entry to the stage is 18729V/unit. 

15 Again the Langmuir-Blodgett relation considered ■ backwards is 
used. The gap spacing that gives a gradient of 18729V/unit 
at the gap entry is p = .007196. The normalised voltage at 
this point is V=0 . 00138794 . To solve for the Decel 2 aft 
face (44) for V=0. To solve for the ground front face (44) 

20 is solved for V=0 . 00138794 • . 

The^ electrodes defined in the previous sections assumed 
negative ions being extracted from 0 up to 27584V. The 
electrodes for the extraction of positive ions are identical 
except that the polarity of the field is reversed. Under 
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this scheme ions are extracted from 27584V down to OV. The 
electrodes are schematically shown in Figure 5. An expanded 
view of the Blocking electrodes is shown in Figure 6. 

DATED this 17th day of October, 2003. 

The Australian National University 

by DAVIES COLL I SON CAVE 

Patent Attorneys for the Applicant 



r,W BUVCT\IJM4JW>pUinj eimikn acvfemriT'lolB 



1 of 4 



V 
+ 



Extractor 

1 



pStagei^ BIock . ng 
I 1 I r- -Stage 2-1 



Stage 1 



Stage 2 



Figure 1 



Extraction 

■ 1 



Stage 1 



Sta'ge2, 




Figure 2 



pweftvvcTunMswpUttM nuxxlcc devfc*nM7/iaro 



2 of 4 




Figure 4 



t .WERWC1MZJM3M ftesnu cxwiion 6wc rtM7/KW)l 



